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Abstract 

We analyze the way in which duality constrains the exact beta function and correlation 
length in single-coupling spin systems. A consistency condition we propose shows very 
concisely the relation between self-dual points and phase transitions, and implies that the 
correlation length must be duality invariant. These ideas are then tested on the 2-d Ising 
model, and used towards finding the exact beta function of the g-state Potts model. Finally, 
a generic procedure is given for identifying a duality symmetry in other single-coupling 
models with a continuous phase transition. 
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Self-duality plays an important role in both statistical mechanics and quantum field theory 
[|TJ . In its simplest manifestation, it relates the partition function of a theory with one coupling 
constant K to the partition function of the same theory evaluated at a different coupling K. 
The most well-known example is that of Kramers- Wannier duality in the 2-dimensional Ising 
model 0, which has various generalizations to different spin models. 

Very recently, self-duality has been a topic of intense studies in the context of field theory 
and string theory. A particularly important case is that of 2-dimensional a-models, where a 
so-called T-duality can be attached to both abelian || and non-abelian [|J isometries of the 
target manifold. In this connection, the focus has most often been on conformally invariant 
backgrounds, because these are the backgrounds of main interest from the point of view of string 
theory. But in these a-model examples duality has been thought to be of wider applicability, 
beyond conformal points. This prompted one of us to investigate the constraints implied by 
duality on renormalization group (RG) flows in these a-models ||. One of the main lessons 
learned from this analysis is that duality forces certain consistency conditions on the a-model 
beta functions. In Ref. f5[ these conditions were shown to independently determine the one- 
loop beta functions (as well as the expected dilaton shift) in a bosonic a-model with an abelian 
isometry. 

In this Brief Report we wish to point out that such consistency conditions are of relevance 
in the context of statistical mechanics as well. For instance, we find that they provide an 
extremely transparent explanation for the close relationship between self-dual points K* and 
phase transitions. They also predict exact duality invariance of the correlation length. In fact, 
starting from a fairly general setting, we find that duality yields a number of results which 
elucidate the RG and phase structure of different models, and this may in turn lead to definite 
predictions on observables in these models. In our effort to test these ideas, we are aided by 
the fact that a full solution of the 2-dimensional Ising model on a square lattice is known [[]], 
while we are on the other hand limited by the scarcity of exact results on most other models 
presenting duality. 

Consider first a theory, for definiteness a spin model, with just one coupling constant K. If 
this theory is self-dual, there exists a relation between this coupling and its dual K, 

k = t(K) , (i) 

such that the partition functions at K and K are related in some simple way, which shall not 
presently be of concern to us. The only essential requirement is that the interaction parts of the 
Hamiltonian and its dual are of precisely the same form, with K replaced by K. We assume 
that the duality map is continuous and 1:1, and that it gives back the original coupling when 
acting twice: 

t(t(K)) = K . (2) 

This implies t'(K) < 0, and we assume K, K > 0. At the self-dual point K* we have t(K*) = K*. 

In quantum field theory, such duality relations can be given meaning only at a given scale of 
the ultraviolet cut-off. In the present spin model context, the ultraviolet cut-off is automatically 
included in the definition of the model (e.g, square lattice, triangular lattice, etc.). We shall 
now explore some simple implications these duality relations have on the RG flow. This flow is 
generated by an infinitesimal change in the ultraviolet cut-off, here as a change in the lattice 
spacing a —>■ a + 5a. In the conventional manner, this leads us to define a beta function by 

dK 

(1{K) = - « or , (3) 
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which stipulates how the coupling constant K is renormalized under a change of scale. 

The consistency condition on (3{K) is arrived at by first applying the operator —ad/ da to 
both sides of the duality relation, Eq. ([[]): 

We now impose that the left hand side of this equation be identified with the beta function itself, 
j3(-), evaluated at K: 

(3{K) = /3(K)t'(K) . (5) 

This is our consistency condition for the beta function. Insofar as K is some function of K, it 
would seem that its RG flow should simply be dictated by that of K. What we are proposing, 
however, is a stronger condition: the requirement in Eq. ([D that the flows of K and K be thus 
"tied together" , as it were, imposes stringent restrictions on what this flow can actually be. A 
simple example below will suffice to show this clearly. 

An immediate result that follows is a general relationship between self-dual points and phase 
transitions. At K = K*, Eq. (|2]) implies that t'(K*) = — 1, and thus there are only two ways of 
satisfying Eq. (|5|) at K = K*\ either 

• the beta function is continuous at K = K*, and (3(K*) =0 , or 

• the beta function is discontinuous at K = K*, with (3(K]K*) = -(3(K [K*) . 

The former case applies to second and higher order phase transitions, or to high/low temperature 
sinks where the correlation length vanishes. The latter possibility, a symmetric discontinuity 
around zero in the beta function, has perhaps a less obvious interpretation. We view it as 
pertaining to first order phase transitions. This is consistent with the discontinuity in the 
correlation length at such transitions, and can also be understood in terms of co-existing phases. 
Namely, we may consider the RG flows on both sides of the first order phase transition, and 
their continuation within one phase only on the other side of the transition. Such analytic 
continuations will typically converge towards a zero close to the phase transition. What duality 
symmetry enforces, then, is that the two continuations occur symmetrically from either side of 
the phase transition, so that the two branches end with equal magnitude and opposite signs just 
at the phase transition point itself. We recall that Kramers and Wannier originally identified 
the self-dual point in the Ising model as a phase transition point by associating it to a point of 
possible non-analyticity in the free energy 0]. We do so in greater generality by associating it 
with a point of onset of scale invariance, with the advantage that our derivation need make no 
assumptions about the lack of phase transition points elsewhere. We now consider what Eq. (|5|) 
entails for the correlation length. 

Under the RG rescaling a — > a' = a + 5a, the correlation length changes according to 
£(K') = t;(K)(l + 5a/ a). Using Eq. (|3|), this gives the well-known relation 

wo - -ffg. (6) 

Together with Eq. (|5|), this then implies that the correlation length should be invariant under 
duality transformations: 

= Z(K) . (7) 
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Although we are interested here mostly in systems with finite order phase transitions, it is clear 
that Eqs. (pj)-([7p can be considered also in the context of systems with more exotic behavior, as 
for instance the 2-d XY model. 

Let us now consider the 2-dimensional Ising model on a square lattice (in the absence of 
external magnetic field) as a specific example. In this case the exact beta function can be 
extracted from Onsager's || original solution (see Ref. 0). In a conventional normalization the 
Hamiltonian is taken to be H = KY^SiSj, where the sum runs over all nearest-neighbor spins. 
The duality relation can then be written 0: 

tanhif = e~ 2k , (8) 

which yields the well-known expression K* = (1/2) ln(l + y/2) = 0.440.... for the self-dual point. 
Moreover, from the exact expression for the infinite-volume correlation length || 

^ = \2K + \ntanhK\ = 2 \K-t(K)\ ' ^ 
one finds, from Eq. (||), the exact beta function: 

, s 1 sinl^i^ , . K-t(K) , , 

As written above, it is easy to see that the beta function does indeed satisfy the consistency 
condition, and that the correlation length is (manifestly) duality invariant. At this point one can 
see more explicitly how restrictive the consistency condition is: using, for instance, Ising duality, 
Eq. (H), if one were to arbitrarily choose some "flow" for K, say K(a) = 1/a (which would yield 
(3{K) = K), then Eq. @, simply expressing a chain rule would, naturally, be satisfied, while 
the consistency condition, Eq. (|), would certainly not. 

When an explicit lattice construction of the duality map is available through an associated 
operator transformation from the lattice to its dual, duality invariance of the correlation length 
can often be checked explicitly from the 2-point correlation function. This is, for example, the 
case of the 2-d Ising model. However, our general considerations are of course not limited to these 
special cases. We find, furthermore, the present derivation to be more direct and transparent. 

While this Ising model solution is completely known, the generalization to the g-state Potts 
model on a square lattice provides interesting ground for conjectures, since there duality is also 
known, but the exact beta function (or correlation length) is not. In these models, the duality 
transformations, K = t q (K), are most simply expressed implicitly through 

(e^-l)(e^-l) = q, (11) 

which reduces to Eq. (||) for q = 2, as it should. In general, for continuous phase transitions with 
a power-law divergence of the correlation length near the critical point K c , we have 

£(K) = A\K-K c \-» + ... , (12) 

and, through use of Eq. (|[) the beta function develops a simple zero, 

P{K) = ^(K-K c ) + ... . (13) 
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This also recovers the well-known fact that the slope of the beta function at K c is determined by 
the critical exponent v = l//3'(K c ). What we have seen is that in the Ising case (where v = 1), 
the behavior above is indeed exact for all couplings if one substitutes \K — K c \ — > \K — K\. If 
we assume that the same also happens in the generalization to q ^ 2, we are led to predict that 
the exact correlation length and beta function are 

UK) = A q \K-t q {K)\^ (14) 

m) = v q l-f q (K) ' (15) 

where v q = (8 — q)/6 for q = 2, 3, 4 . However, one can see that the solution to Eq. (|J), or Eq. (|7|), 
is not unique given only t(K) and the critical behavior at K*, and thus, without further external 
input, we are not quite ready yet to make such a strong conjecture. Furthermore, it is also 
well-known that for q > 5, the transition becomes first order, with a finite correlation length at 
the phase transition. With the beta function as related in Eq. to the correlation length, it 
follows that it must then satisfy the consistency condition with a discontinuity at the critical 
point, which is not the case for f3 q (K) above. We feel, however, that the above considerations 
will be useful in determining exact results, both for g = 3,4 and q > 5, once more information 
becomes available in these models. In particular, in view of these considerations, numerical or 
analytical determinations of these beta functions would be highly interesting. 

We turn now to a question which is motivated by the general character of the Ising model beta 
function and correlation length: as given by Eq. (|9|), C,(K) starts off at zero, grows analytically to 
infinity at K = K*, comes down to the right of that and asymptotes to zero, again analytically. 
Accordingly, the beta function also starts off at zero, becomes negative, crosses zero at K = K*, 
becoming positive and diverging linearly at infinity. We now ask the question: given any other 
model in which the correlation length or beta function behave essentially in this way, is it 
reasonable to expect that it will also present some form of duality? Or, given j3{K) or £(K), can 
one find the duality transformation t(K) uniquely? In fact, one can. We will construct such a 
t(K) here in two alternative ways: first, by a "Taylor reconstruction" procedure, and then by a 
graphical procedure. 

The first procedure consists in simply expanding Eq. (|5|) in a Taylor series around the self- 
dual point. To that end, we first note that whenever £(K) goes to zero or infinity analytically at 
a certain finite K, the beta function has an analytical zero there. This means in particular that 
a Taylor series around the self-dual point makes sense. Moreover, infinitesimally close to the 
self-dual point, K and K are also infinitesimally close. Thus, if we perform such an expansion, 
there will be, order by order in the expansion parameter, a relation giving derivatives of t(K) at 
K* in terms of derivatives of f3(K) and lower derivatives of t(K), also at K*. This procedure is 
straightforward and systematic, and we will not present it explicitly here. The main implication 
is that we can compute the duality map t(K) order by order in a Taylor expansion around K*. 
This (unique) duality map will by construction leave the correlation length £,(K) invariant, while 
the beta function will also satisfy the consistency condition, Eq. ([|). Of course, we cannot from 
this alone argue that t(K) is a genuine duality transformation in the sense of a mapping between 
two dual Hamiltonians (and hence all observables) . But the existence of such a map t(K) in 
this generalized situation is highly suggestive. It would be quite interesting to check by Monte 
Carlo methods whether the resulting map t(K) indeed represents a full duality transformation 
on other observables. One other point is noteworthy: the above procedure cannot be turned 
around to yield (3(K) from t(K), because some derivatives of (3{K) are simply free and drop out 
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of the equations. Of course, this is consistent with the fact that there are many beta functions 
which solve Eq. (0) given t(K), so that there cannot be a unique reconstruction procedure for 

The alternative procedure consists in building t(K) from the £(K) profile: call £,l{K) and 
£r(K) the profiles of £(J£") to the left and to the right of the self-dual point, respectively. If the 
inverse functions £^ and exist, the duality transformation which will respect Eq. fl?]) 
is given by: 

t(K) = &(UK)) , if K>K\ (16) 

= ^(a^) , if k<k\ 

= K* , if K = K* . 

The function t(K) built in this way satisfies all the requirements expected of it and spelled out 
above, and in particular Eqs. (|5]) and (|7]). Since both constructions are unique, they must be 
equivalent to each other. 

There is, of course, still a large number of further applications and tests of the ideas presented 
here. To mention a few, one may for instance extend these considerations to systems which are 
self-dual, but whose dual lies on a different lattice (e.g., 2-d spins on a triangular lattice), or 
yet to systems which are simply not self-dual (e.g., 3-d Ising). It is likely that new features of 
the RG would then come into play, as is for instance the case with scheme dependence in the 
cx-model treated in ||. One may, furthermore, attempt to verify numerically the validity of any 
conjectures resulting from duality considerations. 

Finally, one further context in which duality consistency relations may prove to be extremely 
useful is the quantum Hall effect (QHE). If we assume that the phase structure of the model is 
generated on the (a xx , a xy ) conductivity plane by the modular group SL(2, Z) or an appropriate 
subgroup thereof (cf . |L0f for a detailed treatment of this RG picture) , then the same procedure 
as used above to derive Eq. (|5|) yields immediately some simple but interesting results. For 
instance, the procedure of applying a(d/da) to a duality transformation (in this case a modular 
transformation) and identifying the left hand side with the beta function evaluated at the dual 
point, exactly as above, leads to the requirement that the beta function transform as a modular 
function of weight —2. Such objects are mathematically well-known and studied. Also, if we 
consider that all fixed points are generated by modular transformations of a finite number of 
elliptic points in the fundamental domain, the condition analogous to t'(K*) = —1 (t(K) is 
a modular transformation, K* is an elliptic point) leads to the vanishing of the beta function 
at those points, as expected. It is curious to note that in that case two consecutive duality 
transformations will in general not bring a point back to itself (this is what enforced t'(K*) = 
— 1 for us), and it is instead the existence of elliptic points in the fundamental domain which 
guarantee the vanishing of the beta function at the fixed points. Furthermore, while in the 
single-component case there is a single duality transformation relating two different phases, in 
the QHE, there are an infinite number of duality transformations, relating an infinite number 
of phases. One may hope that such duality considerations, together with particular analyticity 
requirements, could eventually be stringent enough to determine the beta function completely. 

We hope to report on some of these issues in the near future. 

Acknowledgement: We thank M. Asorey and C.B. Lang for discussions. 



5 



References 

[1] J. Kogut, Rev. Mod. Phys. 51 (1979) 659; 
R.D. Savit, Rev. Mod. Phys. 52 (1980) 453. 

[2] H.A. Kramers and G.H. Wannier, Phys. Rev. 60 (1941) 252. 

[3] T.H. Buscher, Phys. Lett. B194 (1987) 59; ibid. B201 (1988) 466. 

[4] X. De La Ossa and F. Quevedo, Nucl. Phys. B403 (1993) 377. 

[5] PE. Haagensen, Phys. Lett. B382 (1996) 356. 

[6] L. Onsager, Phys. Rev. 65 (1944) 117. 

[7] M.N. Barber, in Phase Transitions and Critical Phenomena, Vol. 8 (C. Domb and J.L. 
Lebowitz, Eds.), Academic Press (New York) 1983. 

[8] M.E. Fisher and R.J. Burford, Phys. Rev. 156 (1967) 583. 

[9] F.Y. Wu, Rev. Mod. Phys. 54 (1982) 235; 
H. Giacomini, Phys. Lett. A115 (1986) 13. 

[10] CA. Liitken, Nucl. Phys. B396 (1993) 670, and references therein. Investigations along 
similar lines to the ones sketched here are also being pursued by C. Burgess and C.A. Liitken 
(C. Burgess, private communication). 



6 



